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Abstract 

An exchange energy functional is proposed and tested for obtaining 
a class of excited-state energies using density-functional formalism. 
The functional is the excited-state counterpart of the local-density 
approximation functional for the ground-state. It takes care of the 
state-dependence of the energy functional and leads to highly accurate 
excitation energies. 
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1 Introduction 



Success of density functional theory (DFT) [H |2] for the ground-state calcu- 
lations had prompted search |S1 111 El El IHl 13 UHI for similar theories for 
the excited-states. Over the past decade, time-dependent density-functional 
theory (TDDFT) ^T] has become a standard tool ^21 ^l obtaining 
transition energies and the associated oscillator strengths. However, despite 
it's widespread use, the theory is not without limitations. For example, if 
adiabatic approximations are applied for the exchange- correlation kernel, ex- 
citation energies for double excitation of electrons cannot be obtained J3] 
within TDDFT; obtaining these still remains ^Hl a challenge in the TDDFT 
approach. At the same time, the charm of getting the excitation energy as 
the difference between two total energies remains. This is because one can 
choose the excited-state at will, promoting as many electrons as one wishes 
to a set of chosen orbitals, calculate the corresponding total energy and find 
the excitation energy by subtracting the ground-state energy from it; we 
refer the reader to the works of Ziegler et al , von Barth [H] and Gunnars- 
son and Jones P7| for some of the early attempts to apply this approach to 
study low excitation energies. Thus research in the direction of performing 
a Kohn-Sham like calculation for the excited-states continues. 

A ground-state like DFT approach to obtain the total energy of an excited- 
state has been developed by Gorling [TH] and by Levy and Nagy [12]. The 
theory is based on the constrained-search approach [20] and proposes that 
the energy of an excited-state can also be written as a functional 



of the excited-state density p(r). Here F[p, po] is a bi-density functional 
that depends on the ground-state density po also, and Vexti''!^) is the external 
potential that the electrons are moving in. The bi-density functional for 
the density p of the nth excited-state is defined via the constrained-search 
formulation as 



where \1/ is orthogonal to the lower (n — 1) states of the Hamiltonian, already 
determined by the density po- Such a way of obtaining the functional F[p, po] 
makes it non-universal and also state-dependent. The exchange-correlation 
energy functional -E^.c[p, po] for the excited-state is then obtained by sub- 
tracting from F[p, Po] the non-interacting kinetic energy Ts[p, po] and the 
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Coulomb energy Ucoui[p\ corresponding to p. The non-interacting kinetic 
energy Tg [p, po] is defined in a manner similar to Eq. Q by dropping the 
operator Vee from the right hand side. Thus (for brevity, from here onwards 
we drop po from the argument of the functional) 

E,,[p] = F[p] - Ts[p] - Ucouilp] (3) 

With the assumption that the excited-state density is non-interacting v- 
representable, the density is obtained by solving the excited-state Kohn-Sham 
equation (atomic units are used throughout the paper) 



/ / r — r 



0,(r) = e,0,(r) (4) 



as 

p(r) = S,n,|</.,(r)|2, (5) 

where Uj is the occupation number of the orbital In Eq. the various 
terms have their standard meaning with Vxc{y) representing the exchange- 
correlation potential for the excited-state. It is determined by taking the 
functional derivative of the excited-state exchange-correlation energy func- 
tional. That a Kohn-Sham like calculation can be performed for the excited- 
states was first proposed by Harbola and Sahni j2I] on physical grounds, 
and has been put ^2] on a rigorous mathematical footing recently on the 
basis of differential virial theorem |23] • Calculations of excited-state energies 
based on the Harbola-Sahni work have yielded excellent results fl^ . The 
near exact exchange- correlation potential for the singlet ls2s and triplet 
ls2s excited-state of Helium has also been constructed I^Hl I2Z| recently. 
However, we are not aware of any work where an exchange-correlation func- 
tional for the excited-states has been reported; In performing excited-state 
calculations jUl UH I2H|; either the ground-state functionals or the orbital 
based-theories [2111211 have been employed. The proposition for the construc- 
tion of an excited-state exchange-correlation functional is indeed a difficult 
one since the functional is non- universal and also state-dependent. Thus a 
general functional form for it may not exist. 

Against such a background, we ask if it is at all possible to obtain a 
simple LDA-like functional for the excited-states. To keep matters simple, 
we have been looking at this problem within the exchange-only approxi- 
mation. In this paper we show that it is indeed possible to construct an 
exchange energy functional that gives transition energies comparable to the 
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exact exchange-only theories such as Hartree-Fock [^0], optimized potential 
PU] or the Harbola-Sahni [^1] theory. The construction of the functional is 
based on the homogeneous electron-gas and in finding the final form of the 
functional we are guided mostly by qualitative plausibility arguments. Our 
work is thus exploratory in nature and represents probably the first attempt 
to construct an excited-state exchange-energy functional in terms of the den- 
sity. The evidence of the accuracy of the functional constructed by us is given 
by the results of the transition energies of a large number of excited-states. 
We also refer the reader to ref. |32] for an expression for the change in the 
exchange energy in terms of the ground-state Kohn-Sham orbitals when an 
electron is promoted from a lower energy orbital to a higher one. 

In the present work we take a particular class of excited states in which 
some core orbital are filled, then there are some vacant orbitals and again 
there are some filled orbitals. We construct an LDA-like functional for such 
states in the following section. 



2 Construction of the functional 

As stated above, we now consider such excited-states where the occupation 
of the orbitals is such that the electrons occupy some core orbitals and some 
shell orbitals, leaving the orbitals between the core and the shell regions 
vacant. This is shown schematically in Fig. ([Q). Such an excited-state would 
be obtained, for example, if an electron from the filled orbitals of the ground- 
state is excited to just above the occupied levels. The exact exchange energy 
for a set of occupied orbitals is given as 
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(6) 



So the excited-state exchange energy when an electron is transferred from 
one of the orbitals occupied in the ground-state to the lowest unoccupied 
level is 
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where (prem represents the orbital from which the electron has been removed 
and (padd where the electron is added. The sum over the index j in the sec- 
ond term on the right hand side runs over all the orbitals that are occupied 
in the ground-state and (padd- On the other hand the sum in the fifth term 
runs only over the orbitals occupied in the ground-state. We now attempt 
to make an LDA-like approximation for the excited-state exchange energy 
so that the difference (the last four terms in the equation above) between 
the approximate excited- and ground-state exchange energies is close to that 
given by the exact expression above. In making this approximation accu- 
rate, it is evident that the self-energy terms (third and fourth terms on the 
right hand side of Eq. ((Tj)) for the orbitals (prem and (padd are to be treated 
accurately. 

As the first step towards an excited-state functional, we make the cor- 
respondence between the excited-states that we are considering and similar 
excitations in a homogeneous electron gas (HEG). If the HEG is in it's ground 
state, the electrons are filled up to the Fermi level so that the electrons oc- 
cupy wave- vectors in k — space from k = to kj = (37r^p)3, where p is 
the electron density. On the other hand, in an excited state of the system 
the electrons will occupy k — space differently from the ground state. For 
the kind of excited-states that we consider in this paper, the corresponding 
occupation in the k — space is as follows: The electrons occupy orbitals from 
k = to ki and k2 to k^ with a gap in between as shown in Fig. (^. So the 
excited state density is given by 

p = Pc + Ps (8) 

with 

kl = Sn^Pc (9) 

kl -kl = 37rV™ (10) 

kl-kl = Stt^Ps (11) 

In Eq. (jSI) pc and ps are the core and shell electron density, and in Eq. (fTUj) . 
Pm is the density of the vacant orbitals that lie between the core and the 
shell regions of occupied orbitals. 
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Figure 1: Orbital and the corresponding k — space occupation in the ground 
and the excited state configuration of a homogeneous electron gas(HEG). 



The exchange energy for the HEG that occupies the k— space as described 
above can be obtained exactly and is given as (MLDA stands for modified 
local-density approximation) 



rpMLDA rpcore , rpshell , rpcore— shell 
— ^ + ^X 



where 
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is the exchange energy of the core electrons, 
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is the exchange energy of the electrons in the shell, and 
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(14) 



k2 + kl 
k2 - kl 

(15) 



represents the exchange energy of interaction between the core and the shell 
electrons. Here V is the volume of the HEG. After adding the three terms, 
the exchange-energy can also be written in the form 
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^X 



P [^{h) - ^{h) + e{ki)] dr- J 



.2 ,.2^2 



(^3 ^2 



ks + k2 
ks - k2 



87r3 



{^2 - kiY In 



k2 + kl 
k2 — kl 



dr 



dr 




where e{kf) represents the exchange-energy per particle when the HEG is in 
its ground-state with the Fermi momentum equal to kf. The equation above 
has a nice interpretation: The integral on the right-hand side represents the 
exchange energy of the system of electrons with density p when per electron 
energy is written as [e^k^) — e(/c2) + i-e. the per electron energy is 

given according to the occupation in the k — space (compare with Eq. ()12|) ). 
The log. terms, on the other hand, have no such simple interpretation. They 
have the kinetic energy density in them but we have not been able to write 
the terms in as easy a form as the first term. That the functional above has 
all the right limits if we take ki = k2 or k2 = k^ is easily verified. Finally, 
the modified local-spin density (LSD) functional E^^^^lPa, Pp] in terms of 
the spin densities pa and pp is easily obtained from the functional above as 

i?r"^[p., P,] = \em"'^[2p-] + lE^'^'^m (17) 

Having derived the exchange functional for the HEG, we now apply it 
to the excited-states of various atoms to check if the functional above gives 
exchange energy differences accurately. The excited-states chosen are such 
that they can be represented by a single Slater determinant so that the LDA 
is expected to be a good approximation [31 E] for them. The different radii 
in the k — space {ki, k2 and k^) needed to evaluate the exchange energy 
are found by Eqs. Q, (fTUI) and (fTT|) . For each state (ground and excited), 
the same set of orbitals P2] is employed to get the Hartree-Fock and the 
LSD exchange energies. We calculate the LSD and MLSD exchange energies 
using spherical spin densities since the effect of non-sphericity on the total 
exchange energy should be small [33]. This is because of the fact that in 
the Levy-Nagy formalism ^^1; the excited-state energy is obtained through 
variational minimization. Therefore inclusion of non-sphericity in the density 
would not cause as large a change in the total excited-state energy as the 
use of an appropriate exchange energy functional. Indeed the results for the 
lowest lying multiplets also indicate this |34j . 

In Table 1, we show the difference between the excited-state exchange en- 
ergy and the ground-state exchange energy for some atoms and ions. In the 
first column we give the difference as obtained by the Hartree-Fock expression 
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for the exchange energy. In the second column, the numbers are given for 
both the excited-state and the ground-state exchange energies obtained by 
employing the ground-state LSD functional. The third column gives the ex- 
change energy difference when the excited-state exchange energy is calculated 
using the functional of Eq. (fTTj). It is clearly seen that the ground-state LSD 
approximation underestimates this energy difference. This is not surprising 
since the ground-state functional would give a larger exchange energy for the 
excited-state than what a proper excited-state functional should give. How- 
ever, when the functional of Eq. ()17|) is employed to calculate the exchange 
energy for the excited-states we find, to our surprise, that for the majority 
of the atoms the functional overestimates the differences by a large amount, 
whereas we expected to find the error to be about 10% which is the general 
LDA exchange energy error. We note that this large difference cannot come 
because we have spherical densities. If non-spherical densities are used, the 
difference may increase even further. For example, for the fluorine atom, the 
ground-state exchange energy will become more negative for non-spherical 
densities. On the other hand, the excited-state exchange energy will remain 
unchanged since the density is already spherical. This will result in an even 
larger difference in the exchange energies of the two states. 

We now look for possible sources of error in the exchange-energy differ- 
ences when the functional of Eq. ()17|) is employed to get the exchange energy 
for the excited-states. For this we examine Eq. (|7j) in which the last four 
terms on the right hand side represent the exchange energy difference. Thus 




It is the LDA values to this term that are given in Table 1. The source 
of error in this term we suspect is the LDA treatment of the self-exchange 
energies of the orbitals (prem and (padd involved in the electron transfer. To 
make the functional more accurate we make the self-interaction correction 
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(SIC) for both these orbitals. This is done by subtracting [33] 

= U f Mlimi^rfr,*, + £-"1^(^)1 , (19) 

where p(0) is the orbital density for the orbital 0, from the E^^^^ functional. 
Thus the final expression for the exchange-energy that we have is 

j^MLSDSIC ^ j^MLSD _ _ ^57C[^^^^] (2O) 

This gives the exchange energy difference between the excited-state and the 
ground-state to be 

^Ex = Ex^^^[Pexcited] — Ex^^ [pground] — Ex^[(prem] — Ex^'" [(f)add] (21) 

We have also computed the exchange energy differences given by the func- 
tional in Eq. pO|l and shown them in Table 1. As is evident from the numbers 
displayed there, the functional of Eq. gives highly accurate exchange- 
energy differences for all the systems considered. When the exchange-energy 
difference between the ground- and the excited-state is small, the HP, LSD 
and the functionals derived above, all give roughly the same results. How- 
ever, when this difference is large, the LDA underestimates the magnitude 
of the difference by a large amount whereas the functional of Eq. (jl7p over- 
estimates it. Only when the latter is corrected for the self-interaction then 
the difference is almost the same as the Hartree-Fock difference. Notice that 
SIC is made only for the orbitals involved in the transition. Thus despite 
this explicit orbital dependence, the functional is still quite simple and easy 
to deal with. 

Making the self-interaction correction may deceptively lead the reader to 
momentarily think that our approach may be nothing more than treating 
the exchange energy within the SICLDA approach for both the ground and 
the excited-states. However, this is not so. If the SICLDA exchange energy 
functional is used, the difference between the exchange energies for the two 
states would be (for keeping the expression simple, we are using the same set 
of orbitals for the two states): 

^Ex = E^^^[pexcited] — E^^^[pground] + -E'x^*" [0rem] — E^^^lcpadd] (22) 

Expressions in Eqs. (PT|) and ()22|) differ in two significant ways: First, the 
exchange energy functionals used for the ground and excited states are dif- 
ferent, and secondly E^^'^[(prem] is subtracted in Eq. (PT|) whereas it is added 
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in Eq. (|22j). A careful look at Eq. (j22p also indicates that the excited-state 
energies in SICLDA scheme should not come out to be any different from 
those obtained from the LDA calculations, because E^^^ for the two or- 
bitals involved in the transition would tend to cancel. This is what has been 
observed in the past [TTl ESI EH] • 

We note that we do have a choice of writing the first and the second terms 
in Eq. ^ as 




and then make the self-interaction correction for the orbital (j)rem- How- 
ever, that would bring E^^[(j)j.(.m\ with a positive sign in e^i-sdsic ^ 
the resulting functional will not be as accurate. Although we do not fully 
understand why this happens, we now give a qualitative argument as to why 
the functional of Eq. ()2()|1 gives accurate exchange energy differences. We 
feel that the LDA should be reasonably accurate when the integral over k is 
continuous. As written in Eq. ()18|). the sum in the first term is continuous 
except for the exchange term involving <f>rem and (padd- Thus the LDA to the 
first term should be reasonably accurate. This brings in the self-interaction 
energy of the electron removed with a negative sign in front. By including 
the self-interaction correction for the removed electron only, we find that the 
error in the exchange energy difference reduces to about 10% of the corre- 
sponding HF value. To make the difference even more accurate, we now 
consider the term for the orbital (padd where the electron is added. There, 
when the LDA is made, the electron comes in with its self-interaction so for 
the added orbital too Ex'^' should be subtracted to make the results for the 
energy difference comparable to the Hartree-Fock results. 

Having obtained the functional to calculate accurate exchange energy 
difference, we now apply it to a large number of excited-states of the class 
considered here and find that we get the transition energies very close to 
those given by the Hartree-Fock theory. 
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3 Results 



We employ the exchange functional Ex^^^^^^ proposed above to obtain the 
transition energies for a variety of excitations in different atoms. We find 
that for all the systems the transition energies obtained by us are very close 
to the corresponding Hartree-Fock energies |HZI. Our calculations proceed as 
follows: (a) We get the ground-state energy by solving the Kohn-Sham equa- 
tion with the effective exchange potential calculated using the Dirac formula 
[HHl- (b) We then solve the Kohn-Sham equation with the same (corre- 
sponding to the ground-state) functional for the excited-state configuration. 
This gives us the excited-state energy E^sd, and the LSD exchange energy 
Ejf^ for the excited-state. The difference between E^sd and the ground- 
state energy Eq gives us the transition energy AElsd- (c) We then employ 
the Kohn-Sham orbitals from the excited-state LSD calculation to get the 
modified LSD exchange energy including SIC by employing the functional 
^MLSDSic Qf Eq. (j2nD- (d) The modified transition energy AEf^fXjgDgjQ is 
then given as 

ATP —AT? j_ j^MLSDSIC j^LSD (cyA\ 

i^^MLSDSIC — ^i^LSD + ~ l^^j 

Although we have not performed self-consistent calculations with the new 
functional, self-consistency is not expected to affect the results significantly. 
This is because, as we shall see in the results, the major difference in the 
transition energies given by different functionals arises from the difference in 
the value of the exchange energy itself. We also compare our results with the 
transition energies obtained by the exchange-only time-dependent density- 
funcitional theory (TDDFT) applied within the single-pole approximation. 
We find that our results are comparable to the TDDFT results. In the fol- 
lowing, we have considered three different cases of electron transfer: electron 
making a transition from an 's' to a 'p' orbital; from an 's' to a 'd' orbital 
and from a 'p' to a 'd' orbital. 

3.1 Electron transfer from an 's' to a 'p' orbital 

In this section we consider the cases when one or two electrons are transferred 
from an inner s orbital to an outer p orbital. Shown in Table 2 are the 
transition energies AEhf, ^Elsd and AEmlsdsic ioi some light atoms 
and ions when one of their inner electrons is excited to the lowest available 
orbital. The excitation energy in these systems is such that for some of them 
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^Elsd is close to AEhf but for others it is not. However, AEmlsdsic is 
uniformly accurate for all the systems. We note that the error in AElsd is 
almost fully from the error in the corresponding exchange energy difference. 
This is evident from a comparison of the numbers in Table 1 (for the exchange 
energy differences) and in Table 2. Thus major difference in AE comes from 
the error in calculating the exchange energy. As noted earlier, self- consistency 
effects are much smaller compared to the differences arising from the use of 
the ground-state exchange energy functional for the excited-state also. Our 
results also match well with, and in some cases are better than, the TDDFT 
results shown in the last column of the table. 

In Table 3, we look at the excitation energies of the alkali atoms and Mg'^ 
by exciting an electron from the uppermost orbital to an outer orbital. These 
are weakly bound systems and as such their excitation energies are relatively 
smaller. Thus they provide a good testing ground for the proposed func- 
tional. An interesting point about these systems is that the LSD itself gives 
excitation energies close to the HF excitation energies. It is therefore quite 
gratifying to see that the transition energies obtained by the new functional 
also are of very good quality, although the present method tends to slightly 
overestimate the transition energies. The TDDFT method also gives similar 
numbers although it overestimates the transition energies by a slightly larger 
amount. 

Next we consider some bigger atoms where we can excite the electron 
from more than one inner orbital. Shown in Tables 4 and 5 are the excitation 
energies for the atoms in the third row of the periodic table. In Table 4, we 
consider an electron being excited from the 3s orbital to the 3p orbital. In 
all these case AElsd is smaller than the true energy difference whereas the 
present functional gives highly accurate estimates of the transition energy. 
Notice again that the error in the value of AE^sd arises mainly from the 
error in the exchange energy. The TDDFT results in these cases too are of 
quality comparable to the present method. 

In Table 5, we show the transition energies for the same set of atoms and 
ions as in Table 4, but for the electron now being excited from the 2s orbital 
to the 3p orbital. Consequently the energy of excitation is much larger in 
this case. The LSD in all these cases underestimates the excitation energy, 
whereas the present functional gives accurate results although shghtly over- 
estimating them. However, the error with respect to the LSD is reduced by a 
factor of 5 or more. Thus the proposed functional is accurate for transitions 
from a shallow level as well as from a deep level. We find that the TDDFT 
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results in the present case are not as accurate as in the cases studied above. 

Shown in Table 6 arc the excitation energies for a group of atoms for which 
the LSD gives transition energies very close to the HF excitation energies. In 
all the cases we find that the functional proposed here is able to give accurate 
excitation energies. Thus we find that when the LSD results are accurate, 
so are the results given by the new functional. What is significant, however, 
is that when the LSD results are poor, the new functional properly corrects 
the error in the LSD. 

Finally, we consider the cases where two electrons are excited to the higher 
orbitals. In this case the functional E^^^^^^^ is evaluated by subtracting 
the SIC energy from the E^^^^ for both the electrons. As already pointed 
out, double excitations are difficult to deal with in the TDDFT approach 
to finding excitation energies, because the theory is based on the first-order 
perturbation theory of non-interacting particles. Results for different excita- 
tions for a variety of atomic systems are shown in Table 7. As is evident from 
the table, for all the systems, our method gives excellent results whereas the 
LSD underestimates the energies. In the case of double excitations, no com- 
parison with the TDDFT results can be made because a satisfactory TDDFT 
of double excitations does not exist. 

In all the cases above, we have compared our results with those of Hartree- 
Fock theory and those obtained from exchange-only TDDFT. We do so be- 
cause in our work we have not taken into account the effect of correlations. 
We note that although in atoms Hartree-Fock theory gives total energies 
which are very close to the experimental energies, correlation effects become 
relatively more important in calculating transition energies which is the dif- 
ference between total energies. Thus a comparison with experimental transi- 
tion energies would be meaningful only after correlation effects are properly 
taken into account. 

3.2 Electron transfer from an 's' or a 'p' orbital to a 
'd' orbital 

In this section we consider the case when a 3s or a 3p electron is transferred 
to an incompletely filled 3d orbital in transition metals Sc, Ti, V, Mn, Fe, Co 
and Ni. The results for our calculations on these systems are shown in Tables 
8 and 9. 

In Table 8, numbers are shown for transitions from the 3s orbital of these 
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atoms to their 3d orbital. As is clear from the table, whereas the LSD results 
underestimate the transfer energies, our results compare well with those of 
Hartree-Fock theory. On the other hand, we find that the TDDFT results 
also underestimate the transition energies in comparison to the Hartree-Fock 
theory results. 

In Table 9, results for the transition 3p —>■ 3d are shown. Since there 
are three down spin electrons in the 3p orbital, two electrons are left in it 
after an electron is excited to the 3d orbital. The question arises whether 
to treat these electrons as core electrons, as shell electrons or divide them in 
the core and the shell; the three possibilities give different answers. When 
the electrons are treated as the core electrons, the transition energy comes 
out to be the smallest and it is the largest when they are treated as the shell 
electrons. We have taken the smallest transition energy obtained by us and 
compared it to the energy of transition to the largest possible L value (so 
the smallest transition energy) excited-state for a given spin. In these cases 
too, we find that for Sc, Ti and V, our method gives results which are better 
than the LSD results, but for systems with more than half filled d — shell, our 
method overestimates the transition energy whereas the LSD underestimates 
it. The three cases of Fe, Co and Ni are where our results do not match with 
those of Hartree-Fock theory. 

A tough problem in calculating transfer energies is that of electron trans- 
ferring from the 4s orbital to the 3d orbital in the transition metals considered 
above. The problem has been well investigated f7| EHl in the past and as in 
all the cases considered so far, LSD underestimates these energies by large 
amount. We have applied our functional to obtain these transfer energies to 
see if we could get the correct answer. However, for these s d transfer 
energies, our method gives hardly any improvement over the LSD results; in 
fact for most of the systems, we get a transition energy which is lower than 
the LSD energy. Furher investigations of this problem are being made. 

4 Discussion and concluding remarks 

In the above we have presented a new LDA-like functional for obtaining 
the excitation energies. It has been employed to investigate over 50 excited 
states. The results show that our procedure gives accurate excitation energies 
for all of them, whereas for most of the systems the LSD underestimates 
the energy difference. Thus we have demonstrated that it is possible to 



14 



construct excited-state exchange energy functionals that are capable of giving 
transition energies close to the exact theory. We have worked within the 
exchange-only approximation and have chosen a particular class of excited- 
states. What we have learnt through the study reported here is that there 
exist exchange energy functionals of densities that are more accurate than 
the ground-state functional. However the structure of these functionals is 
not as simple as the ground-state functional. 

In Levy-Nagy theory while defining the bi-density functional for 
the excited-states through constrained search formulation [20J (Eq. (0)), the 
wavefunctions involved in the minimization procedure are those which are 
orthogonal to the lower energy wavefunctions. The latter are supposed to be 
determined by the ground-state density poi^r). In constructing our functional, 
by looking at the orbital occupation in the excited-states, we occupy the 
k — space in a similar manner, representing the unoccupied orbitals by a gap 
in it. Thus in our case the orthogonality condition described above is taken 
care of, to a large extent, by the gap in the k — space. This also reflects an 
implicit dependence on the ground-state density. 

We are now working on functionals for states other than those considered 
in this paper. As pointed out in the introduction, excited-state functionals 
are not universal and therefore have to be dealt with separately for different 
kinds of excited-states. A fundamental question that still remains unan- 
swered is why is it that representing an excited-state by the corresponding 
excited-state of HEG does not by itself give at least as accurate exchange 
energies for the excited-states that the LDA does for the ground-state. We 
do not yet have a satisfactory answer for this. What is clear however is that 
a combination of this functional and appropriate SIC gives highly accurate 
results. In this work, we have also not looked at the correlation energy func- 
tionals. Can correlation energy functionals be developed along similar lines? 
We trust that it should be possible and are working on this problem. 
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Table 1: Difference in the exchange energies of the ground- and excited-states 
of some atoms and ions. The First column gives the atom/ion and the tran- 
sition, the second column the difference AEx^ as obtained in Hartree-Fock 
theory, the third column the difference AEx^^ given by the ground-state 
energy functional. The fourth and the fifth column describes the difference 
as obtained with the functional proposed in this paper. The fourth col- 
umn gives the exchange-energy difference AEx^'^^ obtained by employing 
the functional of Eq. (fT7j) whereas the fifth column gives that given by the 
functional of Eq. ()20j) . AE^^^^^^^ ■ Numbers given are in atomic units. 

atoms/ions AE^^' AE^^'^ A^f^^'^ ^^^mlsusiv 

Li{2s^ '^S 2pi ^P) 0.0278 0.0264 0.0587 0.0282 

B{2s^2p^ 2p ^ 2s^2p2 ^D) 0.0353 0.0319 0.0998 0.0412 

C{2s^2p^ ^ 2s^2p^ ^D) 0.0372 0.0332 0.1188 0.0454 

N{2s'^2p^ "^S 2s^2p'^ *P) 0.0399 0.0353 0.1381 0.0503 

0(25^2/ 3p ^ 2s^2p^ ^P) 0.1582 0.0585 0.2634 0.1624 

F{2s^2p^ 2p ^ 2s^2p^ ^S) 0.3021 0.0891 0.3908 0.2765 

Ne+{2s'^2p^'^P ^2s^2p^^S) 0.3339 0.0722 0.4397 0.3037 

5(35^3/ 3p ^ 3^13^5 3p) 0.1106 0.0475 0.1798 0.1252 

Cl+{3s^3p'^ 3p ^ 3^13^5 3p) 0.1257 0.0483 0.2050 0.1441 

Cl{3s'^3p^^P^3s^3p^^S) 0.2010 0.0603 0.2567 0.1969 
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Table 2: Transition energies, in atomic units, of an electron being excited 
from the 2s orbital of some atoms to their 2p orbital. The first column 
gives this energy as obtained in Hartree-Fock theory. The numbers in the 
second column are obtained by employing the ground-state LDA for both 
the ground- and the excited-state. The last column gives the energies given 
by employing the ground-state LDA for the ground-state and the functional 
of Eq. for the excited-state. 



atoms/ions 

N{2sV 25^2/ 
0+{2s^2p^ "^S 2s^2p^ *P) 
0(25^2/ 3p ^ 2s^2p^ ^P) 
F+(2s22p4 3p ^ 2s^2p^ ^P) 
F{2s'^2p^ 2p _^ 2s^2p^ 
Ne+{2s'^2p^ 2p ^ 2s^2p^ "^S) 



AEhf 


^Elsd 


^Emlsdsic 


AEtddft 


0.4127 


0.3905 


0.4014 


0.4153 


0.5530 


0.5397 


0.5571 


0.5694 


0.6255 


0.5243 


0.6214 


0.5912 


0.7988 


0.6789 


0.8005 


0.7651 


0.8781 


0.6671 


0.8573 


0.7659 


1.0830 


0.8334 


1.0607 


0.9546 



Table 3: The caption is the same as that for Table 2 except that we are now 
considering transitions from the outermost orbital to an upper orbital for 
weakly bound systems. 



atoms/ions AEhf ^Elsd AEmlsdsic AEtddft 

Li{2s^ 2p^ ^P) 0.0677 0.0646 0.0672 0.0724 

Na{3s^ 3p^ ^P) 0.0725 0.0751 0.0753 0.0791 

Mg+{3s^^S ^3p^^P) 0.1578 0.1585 0.1696 0.1734 

K{As^ 4pi 2p) 0.0516 0.0556 0.0580 0.0608 



Table 4: Electron transition energy from the 3s to the 3p orbital in some 
atoms. 



atoms/ions 




AEhf 


AElsd 


AEmlsdsic 


AEtddft 


P(3s23p^ ^5 ^ 
S{3s^3p^ 3p ^ 


3s^3/ 


0.3023 


0.2934 


0.3055 


0.3183 


0.4264 


0.3615 


0.4334 


0.4122 


Cl+i3s^3p^ 3p 


3s^3p^ ^P) 


0.5264 


0.4482 


0.5403 


0.5113 


Cl{3s^3p^ 2p _ 


^ 3si3p6 ^S) 


0.5653 


0.4301 


0.5630 


0.4996 


Ar+{3s^3p^ 2p 


^ 35^3/ ^S) 


0.6769 


0.5174 


0.6766 


0.6007 
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Table 5: Electron transition energy from the 2s to the 3p orbital in the same 
atoms as in Table 4. 



atoms/ions /^Ehf ^Elsd ^Emlsdsic AEtddft 

6.4188 6.9564 6.1573 

7.7337 8.3271 7.4533 

9.2551 9.8997 8.9618 

9.1653 9.8171 8.8686 

10.8009 11.5061 10.1901 



Table 6: Electron transition energy when the upper state is not the lowest 
energy multiplet. 



atoms/ions 




AEhf 


AElsd 


AEmlsdsic 


AEtddft 




■ 2s'2p^ '^D) 
25^2/ ^D) 


0.2172 


0.1993 


0.2061 


0.2168 


C(2s22p2 3p ^ 


0.3290 


0.3078 


0.3216 


0.3325 


■ 2s^2p^ ^D) 


0.2942 


0.2878 


0.2967 


0.3090 


N+{2s^2p^ '^P - 


2s^2p^ '^D) 
3si3p2 -^D) 


0.4140 


0.4149 


0.4305 


0.4433 


Si+{3s^3p^ 2p 
Si{3s''3p'' 3p - 


0.2743 


0.2632 


0.2799 


0.2864 


> 3s^3p^ ^D) 


0.2343 


0.2356 


0.2442 


0.2567 



P{2s^3p^ 2s%/ *P) 6.8820 

S{2s'^3p'^ 3p _^ 2s^3p^ 3pj 8 2456 

C/+(2s23p^ 3p _^ 2s^3p^ 9.8117 

Cl{2s^3p^ 2P ^ 25^3^6 ^S) 9.7143 

^/•+(2.s23/ -P 2.s^3/ -S) 11.3926 
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Table 7: Excitation energies of some atoms when two electrons are excited. 



atoms/ions 

Be{2s' 2f ^D) 

C+{2s^2p^ 2p ^ ^D) 

7V+(2s22]92 3p ^ 2p4 

C'+(2s22p3 4^ ^ 2/ 2p) 
0(25^2/ 3p ^ 2/ i^S) 
F+(2s22p4 3p ^ 2p6 ^5) 

,S(3s23/ 3p ^ 

P(3s23p3 4^ ^ 3p5 2p) 
52+(3s23pl 2p ^ 3p3 2p,) 

5^(3s23p2 3p ^ 3p) 



AEhf ^Elsd ^Emlsdsic 

0.2718 0.2538 0.2655 

0.4698 0.4117 0.4798 

0.6966 0.6211 0.7180 

0.7427 0.5950 0.7312 

1.0234 0.8369 1.0143 

1.1789 0.9440 1.1785 

1.5444 1.2552 1.5480 

1.5032 1.1333 1.4736 

1.8983 1.4381 1.8494 

0.2578 0.2555 0.2651 

1.0273 0.7807 1.0266 

0.8539 0.6927 0.8680 

0.5856 0.5377 0.6230 

0.5860 0.4928 0.5986 

1.2535 0.9551 1.2516 



0.2718 
0.4698 
0.6966 
0.7427 
1.0234 
1.1789 
1.5444 
1.5032 
1.8983 
0.2578 
1.0273 
0.8539 
0.5856 
0.5860 
1.2535 



Table 8: Electron transition energy when an 's' electron is transferred to a 
'd' orbital. 

atoms/ions 

Sc{3s'3d^ 3s^3d'^ ''G) 

Ti{3s^3(f 3p ^ 3si3rf3 5p) 
Ti{3s''3<P ^ 3s^3d^ ^H) 

y(3s23rf3 4p ^ 3sl3(/4 4^) 

Mn{3s'^3d^ 3s^3rf6 ^D) 

Fe{3s^3d^ 3s^3d.^ ^F) 

Co{3s^3d'' ^F 3s^3d^ ^F) 

Ni{3s^3d^ ^F 3s^3d^ ^D) 



AEhf 


AElsd 


AEmlsdsic 


AEtddft 


2.1562 


1.8584 


2.1223 


1.8649 


2.2453 


1.9740 


2.2061 




2.3861 


2.0827 


2.3649 


2.0951 


2.6098 


2.3107 


2.6106 


2.3266 


3.1331 


2.7860 


3.1199 


2.8062 


3.4187 


3.0483 


3.4527 


3.0755 


3.7623 


3.3178 


3.7955 


3.3516 


4.1204 


3.5949 


4.1476 


3.6351 
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Table 9: Electron transition energy when a 'p' electron is transferred to a 'd' 
orbital. 



atoms/ions ^Ehf ^Elsd ^Emlsdsic AEtddft 

Sc{3p^3d^^D ^3p^3(P^H) 1.1295 1.1018 1.1245 1.2128 

^^(3/3rf2 3f _^ 3p^3d^ ^I) 1.2698 1.2478 1.2728 1.3586 

V{3p^3d^ *F 3j9^3d^ ^/) 1.4153 1.3959 1.4227 1.5042 

Mn{3p^3d^^S ^3p^3d^^F) 1.7270 1.6431 1.6726 1.8073 

Fe(3/3(^6 ^ 3p^3(i^ ^G) 1.8785 1.8784 2.0061 1.9898 

Co(3/3d^ ^ 3p-^3d^ ^G) 2.1178 2.0568 2.2778 2.1755 

Ni{3p^3d^^F ^3p^3d^^F) 2.4232 2.2402 2.5518 2.3656 
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